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^ : ABSTRACT 

O I We present a C/tan(ira observation of the H II region RCW89. The nebula lies on 10' 

^ I north from the central pulsar PSR B1509-58, and it has been suggested that the nebula 

' is irradiated by the pulsar jet. We performed a spectral analysis of the seven brightest 

■ emitting regions aligned in a "horse-shoe" shape, and found that the temperature of the 
^ ' knots increases along the "horse-shoe" in the clockwise direction, while, in contrast, 

■ the ionization parameter n^t decreases. This strongly supports a picture of energy 
transfer via the precessing pulsar jet. We examined the energy budget assuming that 
RCW89 is powered by the pulsar jet, and confirmed that the pulsar rotational energy 
loss is sufficient to drive the nebula. The rate of energy injection into RCW89 by 
the jet was estimated from the synchrotron radiation flux. We obtained a heating 
time-scale of 1400 yr, which is consistent with the pulsar characteristic age of 1700 yr. 
To explain the temperature gradient, we discuss the cooling process for plasma clouds 
in RCW89. We argue that the plasma clouds can be cooled down by the adiabatic 
expansion within 70 yr, and form the temperature gradient reflecting the sequential 
heating by the precessing pulsar jet. We also determined the velocities of the individual 
plasma clouds by spectral fltting. The plasma clouds in RCW89 are moving away at 
240 ~ 860 km s~^, which constrains the inclination angle of the pulsar spin axis i > 50° 
and the expanding velocity of the shell as f shell > HOO km s~^. 



Subject headings: pulsars, supernovae, supernova remnants 
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1. Introduction 

Recent high-resolution X-ray imaging of young energetic pulsars have revealed the ubiquitous 
presence of "jets" or collimated outflows along the symmetry axis, as seen in the Crab pulsar 
(Brinkmann, Aschenbach, and Langmeier, 1985; Hester et al. 1995; Weisskopf et al. 2000), the 
Vela pulsar (Helfand, Gotthelf, & Halpern 2001; Pavlov et al. 2002), and PSR B1509-58 (Trussoni 
et al. 1996; Gaensler et al. 2002). These jets are probably formed by the relativistic outflows of 
pulsar magnetospheres, but their composition, energy contents, coUimation mechanism, and their 
relation to the ring- or disk-like pulsar wind nebulae have not been understood. Investigations of 
pulsar jets allow us to understand the mechanism of particle acceleration and interaction between 
the pulsar wind and the interstellar medium. Particularly intriguing is the case with PSR B1509- 
58, which is accompanied by an unidentified thermal nebula RCW89 on the terminus of the pulsar 
jet. Because of its location, RCW89 has been considered to be a contact point of the jet and the 
surrounding SNR or the ambient medium. If it is indeed the case RCW89 could be a precious 
probe for the mysterious pulsar jet. 

The supernova remnant (SNR) G320.4-01.2 (MSH15-52, Kes23) has bilateral radio shells with 
a diameter of ~ 30 arcmin (Caswell, Milne & Wellington 1981). PSR B1509-58, a 150 ms X-ray, 
radio, and gamma-ray pulsar, is located at the center of G320.4-01.2 (Seward & Harnden 1982; 
Manchester, Tuohy, & D'Amico 1982). Prom the spin parameters, a characteristic age Tc = 1700 
yr (Kaspi et al. 1993), a spin-down luminosity £■ = 1.8 x 10^^ erg s~^, and a surface magnetic 
field Bp = 1.5 X 10^^ G have been obtained, making it one of the youngest, the most energetic, 
and the highest field pulsars ever known. The age of 6—20 kyr for the SNR has been derived 
by assuming standard parameters for the ISM and for the supernova explosion. But this age is 
an order of magnitude larger than the characteristic age of the pulsar. Although the SNR has a 
complex structure with multiple components, H I absorption measurements confirm that it is a 
single SNR at a distance of 5.2±1.4 kpc (Gaensler et al. 1999), suggesting that both the pulsar 
and the SNR originate from the same supernova. 

The Ha emission nebula RCW89 was discovered by Rogers, Campbell, & Whiteoak (1960) 
and coincides with the brightest portion of G320.4-01.2 in radio. It also coincides with one (north 
nebula) of the two clearly distinguishable X-ray nebulae observed by the Einstein Observatory 
(Seward et al. 1983). The X-ray images taken by RO SAT show the distribution of the plasma 
clumps into a "horse-shoe" like shape (Trussoni et al. 1996; Brazier et al. 1997). Two possible 
candidates have been proposed for the energy source of the thermal nebula: the supernova blast 
wave, and the pulsar jet (Manchester 1987). Prom the analysis of ASCA data, Tamura et al. (1996) 
found a bridge of non-thermal emission connecting the pulsar to the thermal emission. They also 
found that in RCW89 region thermal emission with prominent line features is dominant. The 
analysis of the thermal X-ray spectra showed that its energy content was consistent with being 
powered by the pulsar through the jet. They also argued that a less collimated jet or the precession 
of the jet is required for heating the spatially extended plasma. 

Gaensler et al. (1999) reported detailed radio observations using the Australia Telescope 
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Compact Array (ATCA). Due to the correspondence of the radio synchrotron knots with the 
clumps in the X-ray image, they propose that the pulsar is interacting with RCW89 through the 
jet, and that the non-thermal knots are embedded in the diffuse nebula emitting the thermal 
X-rays. This picture was however questioned by the Chandra data showing that all the knots in 
RCW89 have thermal spectra (Gaensler et al.2002). Thus the model involving thermal clumps 
embedded in a diffuse synchrotron nebula (as proposed by Tamura et al. 1996) seems more likely. 
In this case the energy budget discussed by Tamura et al. (1996) must be reconsidered; Chandra 
has revealed that the energy content in the knots is by two orders of magnitude smaller than that 
deduced from the ASCA observations. 



2. Observation and Results 

G320.4-01.2 was observed on 2000 August 14 with the Advanced CCD Imaging Spectrometer 
(ACIS) I-array aboard Chandra. The net exposure time was 19039 s after removing small intervals 
in which no data were recorded (Gaensler et al. 2002). A smoothed X-ray image in the energy 
range 0.4—8.0 keV is shown in Figure 1. 

PSR B1509-58, accompanied by the bright PWN, is at the center of the FOV of the ACIS-I 
and the jet features are extended toward the north-west and south-east direction. RCW89 hes 
on the end point of the north-west jet, exhibiting a complicated clumpy structure. We extracted 
photons from the brightest seven emitting regions (labeled as 1—7 in clockwise direction in Fig. 
1) and the jet region for comparison, using the psextract script of CIAO 2.3 and carried out a 
spectral analysis. The background was taken from the area marked in Figure 1. Figure 2 shows 
the spectra of the north jet and region A. The spectrum of the north jet is well fitted by a single 
power law function with a photon index of F = 2.07 ±0.11 and a hydrogen absorption column 
density of Nr = (0.86 ± 0.09) x 10^^ cm-^ with a x^red = 0-981 70 d.o.f.. The absorption 
corrected flux density for the energy range 0.5—5.0 keV is (1.26 ± 0.09) x 10~^^ ergs s~^. 

The energy spectra of RCW89 regions are quite complex and show strong emission lines as 
well as a hard high energy tail, which has already been found in the ASCA data (Tamura et al 
1996). For the analysis of individual knots in RCW89, the photon statistics was not sufficient 
to determine the column density Nh and the photon index F for each knot independently. To 
determine the spectral parameters for each knot we first analyzed the spectrum of the total region 
A in Figure 1, which contains all of the bright knots 2—7. The spectrum of the region "A" is 
shown in Figure 2. Wc applied two component model consisting of a non-equilibrium ionization 
model ("VNEI" in XSPEC 11.3.0) for the thermal component and a power law for the high energy 
tail. In the NEI model, the metal abundances of Ne, Mg, Si, and Fe, which all have remarkable 
line emission, were set to be variable. The best fit parameters arc shown in Table 1. The column 
density and the photon index were obtained as Nh = (1.18±0.01) x 10^^ cm~^ and F = 2.48±0.05 
with a Xred ~ 1-727 for 227 d.o.f. The fit is rather poor which is mainly due to large residuals 
in the energy range of 1.0—1.3 keV. For the individual spectral analysis we used fixed values of 
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Nh — 1.18 X lO^^cm ^ and F = 2.48. The results of the spectral analysis for each knot are shown 
in Table 2. 

Figure 3 shows the temperatures and ionization parameters ngt from the individual fits to 
the knots of RCW89. This allows a crude estimate of the physical age of each cloud. The 
temperatures tend to increase along the "horse-shoe" in a clockwise direction from knots 1 to 7, 
while the ionization parameter decreases. This result implies that the knots on the north edge of 
RCW89 consist of freshly heated hot plasma and are still in the state of non-equilibrium ionization 
while the clouds in the southern regions are close to ionization equilibrium. 

3. Discussion 

In the previous sections, we show clear evidence for variations of plasma parameters (temper- 
ature and age) in bright knots in RCW89. There arc two possibilities to account for the thermal 
X-ray emission from RCW89, (1) shock heating due to the supernova blast wave and (2) heating 
by the pulsar jet of PSR B1509-58. To constrain the energy source of RCW89, we evaluate the 
ages of the plasma clouds. The Chandra image yields the volumes for emitting plasma clouds, 
assuming that the shapes of the clouds are ellipsoids with major and minor axes of la and in 
arcsec. We obtain the number densities of the electrons in the clouds as 

where EM is the emission measure for each cloud derived from the spectral fits above. Dividing 
the ionization parameter Ugt by the electron number density Ue provides the ionization age t-p. 

The resulting ages of the plasma clouds arc shown in Table 3. Their ages are in the ranged of 5 
~ 700 yr, which is clearly different from each other. Thus, it is not likely that all the knots of 
RCW89 were heated simultaneously by the supernova blast-wave. 

Next, we examine the second hypothesis (that the knots were heated by the pulsar jet) based 
on the energy budget of the plasma clouds. The internal thermal energy stored in a plasma cloud 
is estimated as 

E., = 6.7 X 10- (/. ( (—) (^Y" ergs, (3) 

which are shown as Eth in Table 3. The total energy in the seven knots amounts to 2.6 x 
10^^ (o?L/5.2kpc)^^^ ergs. The heating time scale can be constrained by taking the observable pulsar 
jet flow into account. The non-thermal flux in the north jet region represents the synchrotron 
emission from the relativistic electrons. From the spectral fit in Table 1, we therefore infer that 
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unabsorbed spectral luminosity at 1 keV is Li ]jev = 2.1 x 10^^ ergs s~^ keV~^ for a distance of 
= 5.2 kpc. Wc assume the jet to be a cylinder of length Ij = 3 x 10^^ cm and cross section 
Sj = 2.0 X 10^^ cm^. Assuming equipartition between electrons and magnetic field, we obtain the 
minimum energy of electrons and magnetic field in the jet, 

Ej^t = 6.4 X 10^^(1 + iif'^fl'' (cos i)'^!^ ergs, (4) 

where / is the filling factor, and /i is the ratio of ion to electron energy, i is the inclination angle 
of the jet. The obtained energy E,j(.t corresponds to a magnetic field Bj^.^ ~ 34/iG in the north 
jet, which is consistent with B ~ 25/iG in the south jet evaluated by Gaensler et al. (2002). The 
observed photon index for the north jet is ~2.0, while the south jet in the vicinity of the pulsar has 
a photon index of 1.6. This result implies that the synchrotron electrons emitting in the observed 
energy band had been cooled before the jet reached the terminal. In a magnetic field B ~ 34/iG, 
the life time of synchrotron electrons which emit 1 keV photons is ~ 50 yr. Assuming that the 
electrons are accelerated in the vicinity of the pulsar, the life time of the electron constrain the 
velocity of the north jet, 

where / is the distance between the "jet (north)" region marked in Figure 1 and the central pulsar. 

For a bulk velocity of the jet f blk ~ 0.5c, analogous to that of the Crab nebula(Hester et al. 
2002), the power injected into RCW89 amounts to 

= 3.2 X 10=''^^(1 + ^I'fl' (cos^)'^/^ (^) ergs s-^ (6) 

This corresponds to 1.7% of the pulsar's total spin-down luminosity. If the ratio at which the 
pulsar supplies energy to the north jet had not changed, we can evaluate the initial rotation 
period as 35 ms for a pulsar's moment of inertia 10'^^ g cm^. With a constant braking index of 
n = 2.8, a pulse period of P = 150 msec, and its first derivative P = 1.5 x 10~^^ s s~^ (Kaspi et 
al. 1993), the required time to heat RCW89 is then 1400 yr. The ages of oldest plasma estimated 
from the ionization time scale and the heating time scale required by the energy budget are both 
of the order of ~ 10^ yr, in a good agreement with the pulsar's characteristic age of 1700 yr. Thus 
our results favors the scenario of heating by the pulsar jet rather than the simultaneous shock 
heating by the blast wave. 

With the obtained magnetic field, we can also estimate the luminosity in TeV range. In the 
magnetic field of 34 /xG, the Lorentz factor of the electrons which emit synchrotron photons at 
1 keV is ~ 2.3 x 10^. Such ultra-relativistic electrons can up scatter CMB photons up to ~ 1 
TeV gamma-ray. The ratio of the inverse Compton luminosity to the synchrotron luminosity can 
be described as Lie /L sync ~ ucMs/umag, where ucmb and Umag are the energy densities of the 
cosmic microwave background and that of the magnetic field, respectively. Here, we must compare 
the synchrotron/inverse Compton luminosity correspond to the same population of relativistic 
electrons. The observed synchrotron luminosity in 0.5-5.0 keV is Lgync 4.0 X 10 ergs s . ^Ve 
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infer a luminosity of the inverse Compton component in the north jet Ljc ~ 3.0 x 10^^ ergs s""^ 
in 1.3-13 TeV which corresponds to the energy of CMB photons scattered by the synchrotron 
electrons emitting in 0.5-5.0 keV. Since the expected energy flux of 1.0 x 10~^^ ergs s~^ cm~^ is 
lower than the detection limit of HESS, we cannot expect to detect the TeV gamma-rays from the 
north jet. 

The observed temperature of the elder plasma clouds are cooler. This result implies that 
an effective cooling process was working and have formed the observed temperature gradient. 
To investigate how the temperature gradient was formed we consider the cooling process of the 
plasma clouds. 

First we examine the radiative cooling. The dominant emission processes from RCW89 are 
bremsstrahlung and line emission. However, the radiation fluxes observed by Chandra are too 
small to explain the temperature gradient. We find that the cooling time scale of the knots due 
to optically thin thermal radiation is longer than ~ 4000 yr which is much longer than the pulsar 
age, and cannot explain the temperature gradient. 

Another plausible process is an adiabatic cooling. We discuss two possible cases of adiabatic 
expansion. In the first case {^^pure expansion"), the heated gas volume expands without mixing 
with the ambient matter. In the second case {''SedoiJ^), the hot gas expands like in the Sedov 
phase of an SNR, where the expanding shock front mixes and heats the ambient matter. 

In the "]9ure" adiabatic expansion, the temperature T and the volume V follow the relation 
j^yj-i — const, (7 is the specific heat ratio; 7 = 5/3 for a monoatomic gas). We consider, for 
example, the case of the plasma cloud labeled "7" in Figure 1 with a radius of ro = 2 x lO^'' cm 
and a temperature of To = 0.6 keV. From the V — T relation, a radius of r = 2.8 ~ 10^^ cm is 
expected when the knot cools down to a temperature of 0.3 keV. Assuming that the expansion 

/ N 1/2 

velocity equals to the sound velocity in the plasma, = ( ^ 1 and the thermal pressure of 
p — 2nekT, the expansion velocity can be expressed as a function of radius r as 

dr _ f 2^kTorl V'^ 1 
dt \ nip J r 

Integration yields the time required to expand from ro to r as 

From this equation, the cooling time is estimated as tcooi ~ 70 yr. In this calculation we have 
assumed that the density of the plasma cloud is much higher than that of the ambient ISM, which 
thus does not affect the expansion velocity. 

Next, we examine the "^erfoz/' case. The radio observations showed that RCW89 lies on the 
north shell of the SNR (Gaensler et al. 1999, 2002). Therefore, the jet may heat a small volume 
of the north shell to form a hot plasma which expands into the ambient matter. In the Sedov 
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model the injected energy Eq is divided into internal energy and into kinetic energy as 



(9) 



where p is the pressure, p is mass density, and v is the expansion velocity. We crudely approximate 
the velocity as v ~ c^, which yields 



15 \ /Eo} 



(10) 



where Eq is the injected energy. For a radius of r = 2 x 10^^ cm, velocity was estimated to be 
V ~ 400 km s~^. The radius r can be described as a function of time t by integrating the velocity 



r = 1.0('^V^^^/^ 



The relation between the temperature kT and the radius r is obtained by substituting the equation 
of the state p — 2nekT into the Eq. (10) as 
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(12) 



where k is Boltzmann constant and Up and are the number densities of proton and electron, 
respectively. Then we obtain the time scale to cool from kTg to kT as 



t. 



cool 
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5/6' 



(13) 



For knot- 7 {tq = 
estimated as Eq 
40 years. 



2 X 10^'' cm, T ~ 0.6 keV, rie = 130cm ) the initially injected energy was 
1.3 X 10"^^ ergs. The time required to cool from 0.6 keV to 0.3 keV is about ~ 



Both the "pure expansion" and "Sedov model" can provide sufficient cooling to form the 
temperature gradient within a heating time scale. At this moment it is hard to say which process 
is more appropriate since we do not understand well enough the environment surrounding the 
plasma clouds. Moreover, the knots are embedded in a diffuse synchrotron emission. It may be 
therefore possible that these plasma clouds are confined by the ambient magnetic field, and their 
expansion may not be as free as we have assumed here in the discussion. 

The differences in the temperatures and plasma ages of the knots indicate that RCW89 was 
heated in sequence. If these plasma clouds have been heated by the pulsar jet, the direction of the 
jet would vary with the period. Prom the tendency of the temperature and the plasma age, we 
infer that the interaction point between the jet and the ISM has shifted along the "horse-shoe" 
in clockwise direction. This implies the "precession like" motion of the jet. Preceding discussions 
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about the energy budget and the coohng process indicate that the jet has made only one precession 
since the birth of the pulsar within a time scale of 1000 yr. 

Considering the origin of the matter which forms the hot clouds, the abundance of elements 
may provide a precious clue. Prom the spectral fitting we have obtained the metal abundances of 
the knots in RCW89. The observed iron abundances tend to be smaller than the solar abundance, 
while Ne and Mg are relatively rich as shown in Table 2. This suggests that the plasma in RCW89 
may originate from a type-II supernova. If RCW89 consists of the ejecta from the supernova 
associated with the central pulsar, the ejecta must have traveled at least ~ 7.5 pc from the SNR 
center. In order to arrive at RCW89 within 1700 yr, the required velocity of the ejecta is more 
than ~ 4000 km s~^, which is rather fast for the standard ISM condition. 

The fits to the thermal emission component of the knots required non-zero velocity of the 
emitting gas to account for slight shifts of the fitted line centroids. We can constrain the geometry 
and the velocity of the ejecta (RCW89) using the velocity of the plasma clouds. Figure 4 (left) 
shows the observed redshift for the each knot. Except for the knot 3 which has the poorest photon 
statistic, the redshifts he on a smooth curve from 0.9 x 10~^ to 2.9 x 10~^, which correspond to 
the velocities of 240 ~ 860 km s~^ along the line of sight. Now we can consider the geometry as 
shown in Figure 4 (right). The SNR system is approaching us at ~ 70 km s~^ due to the rotation 
of the Galaxy. Using the inclination angle i between the pulsar spin axis and the line of sight and 
the half cone angle a of the precession of the jet, the apparent velocity Vapp is described as 

^^app max = ^^shell ■ cos(i - a) - 70° ( km s~^), (14) 
^^app min = t^shell • COs(i + a) - 70° ( km s~^), (15) 

where t'shell is the velocity of the shell boarding RCW89. Dividing eq. 14 with eq. 15 to cancel 
■^shell yields the relation of i and a as. 

cos(/ - a) _ t;app max + 70° 
COs(i + a) Vg^pp + 70° 

The observed X-ray image (Figure 1) indicates that the open angle a is smaller than 20°. Thus 
the inclination angle is constrained as i > 50°, which is consistent with the estimation i > 70° by 
Brazier & Becker. (1997). If we know i and a, we can calculate the velocity of the shell. Using 
equation 14 and 15 give 

Vapp max — "^app min ,^ „n 

^^l^^ll ^ 2sinisina " ^^^^ 

For the inclination angle i > 50° (a < 20°), the expanding velocity of the shell is obtained as 
■^shell > 1100 km s~^ (For f shell ~ 4000 km s~^, i ~ 80° and a ~ 5° are required ). 



4. Summctry 



We presented a Chandra study of the H II region RCW89. The observed X-ray image (Figure 
1) revealed the complicated clump structure distributed along a "horse-shoe" shape. To investigate 
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each plasma cloud, we selected the 7 brightest knots marked in Figure 1 and performed spectral 
analysis for each knot individually. In the model fitting we assumed a NEI model and power law 
function. 

The resulting temperature tend to increases along the "horse-shoe" in the clockwise direction, 
while the ionization time scale r decreases. The gradation of temperature and r indicate that 
the knots on the north edge of RCW89 are freshly heated hot plasma and are still in the state 
of non-equilibrium ionization, while the plasma in the south region are older and more close to 
ionization equilibrium. This new finding is consistent with the picture that RCW89 nebula is 
heated by the pulsar jet sequentially. 

To confirm the hypotheses, we examined the energy budget for RCW89. The observed emis- 
sion measure and volume yielded stored thermal energy in the knots as Eth = 2.6 x 10'^'' ergs. On 
the other hand, the pulsar jet is supplying relativistic electrons and magnetic energy of 3.2 x 10"^^ 
ergs s~^ into RCW89 assuming a condition of equipartition and a bulk velocity of the jet of 
Vbi.k — 0.5c. The jet requires ~ 1400 yr to heat the knots, which is consistent with the pulsar 
characteristic age. 

Next we discussed the cooling process to explain the gradation of the temperature. The 
most plausible process is adiabatic expansion. If the expansion velocity is the sonic velocity in 
the thermal plasma, it takes ~ 70 yr for cooling down the knot-7 with typical values of internal 
energy Eth = 1.3 x 10^^ ergs and radius of r = 2 x 10^^ cm from 0.6 keV to 0.3 keV. It is sufficient 
to form the temperature gradient within the pulsar time scale of 1700 yr. 

The spectral fitting provides information on the metal abundance of each knot. The Fe- 
poor spectra indicate that the matter consisting RCW89 originated in a Type-II supernova. It is 
consistent with the existence of the pulsar. 

The obtained redshifts of the each knots vary with their locations. We estimated the geometry 
of the pulsar and the velocity of the shell of the SNR using the Doppler shifts of the spectra and the 
X-ray image. With the precession half cone angle of a < 20°, the inclination angle is constrained 
as i > 50°, and the velocity of the shell as I'shell > l^OO km s~^. 
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Table 1. Spectral fits to region 


"A" and "jets' 


' in Figure 1. 




Parameters 


region-A 


jet (north) 


jet (south) 


Column density (lO^^cm"^) 


1.18 ±0.01 


0.86 ± 0.09 


1.02 ± 0.08 


Thermal component 

Temperature (keV) 

ngt'^ (10^^ s cm"^) 
EM^ {10^^ cm-^) 


0.38 ±0.01 

0.71 ±0.05 
3.42 ±0.03 






Abundance (solar abundance) 
Ne 
Mg 
Si 
Fe 


1.00 ±0.03 
0.92 ±0.02 

0.41 ±0.04 
0.32 ± 0.02 






Non-thermal component 
Photon index 

Fluxo.5-5.Okev (10~^^ ergs s~^ cm-^) 


2.48 ±0.05 
8.67 ±0.05 


2.07 ±0.11 
1.26 ±0.09 


1.66 ±0.09 
4.16 ±0.26 


xl (^) 


1.727 (227) 


0.981 (70) 


0.493 (172) 



Note. — Errors are with 1 a confidence. 

^Ionization timescale of the non-equilibrium ionization plasma. 

"^Emission measure J fienndV {di/ 5.2^.^0)"^ ( is the distance to RCW89). 



Table 2. Results from spectra fitting 



Region ^ kT 


net EM F^Yi 






Abundance 




Rcdshift 


xl {y) 


(keV) 


(10"cm-3s) (1056pin-3) (xlO"") 


(xlO-") 


Ne 


Mg Si 


Fe 


(xlO-3) 





1. 

2. 
3. 
4., 
5. 
6. 
7. 



21+" °^ 

"•^"-0.01 



0.20^ 



F0.04 
-0.03 

41+0 02 

"•^-^-0.01 

-,+0.01 
-0.01 
+0.02 



o.4o: 

0.52 



-0.02 

n 67+0 08 

^•"'-0.08 



> 11(24.2) 

> 11(28.8) 
1.46+1-59 



-0.68 
1-0.16 
-0.13 

70+0 07 

+0.06 
-0.05 
+0.08 
0.05 



1.07 



0.46 
0.23 



3-94_2.28 
9 '?n+0-64 
^•'-'"-0.52 

11.561^^^ 

h0.16 
-0.04 
fO.28 
-0.29 
7+0.06 
-0.12 

t^s+o-i^ 



1.43] 
8.22^ 
1.17"' 



1 03+114 

1 2q+o-''i 

-^•^^-0.20 
-^^•^"-0.32 



+0.53 
0.48 
+0.48 



9.89 
2 94 

^■''^-0.38 



9 OO+0.30 
^••^•^-0.34 
4 17+0.42 
^•^ ' -0.42 



0.41 
0.41 
+0.34 
0.42 
+0.59 
-0.59 
+0.47 



5.33 
2.01 
3.11 

2-96_o.46 

o 4^+0.64 
^•^"3-0.80 



1 64+0-24 

l^-D^-0.24 
1 94+0.22 
^•^^-0.23 
1+0.14 
-0.15 
1 (-9+0.13 
-•-■^^-0.14 

1 -ro+0.05 

-L--L«-o.05 

F0.15 
-0.14 

2 16+0-48 

^•-'-'-'-0.29 



0.54^ 



2.16^ 



2 46+0-68 

^■+"-0.36 

2 qq+0.61 

^•^^-0.30 

1 02+0-87 

-^•"^-0.19 

1 46+0-13 

^■^"-0.11 

0.93t"°^ 



-0.04 
hO.lO 
-0.11 

1 i^'^+o-ii 

l.JO_g 23 



1.77^ 



1 02+104 

l-"^-1.02 

< 0.90(0.40) 



1 16+0-53 

^•-^"-0.56 
;+0.13 
-0.12 

fSO+0-06 

1+0.11 
-0.11 

0.38t°;?i 



0.35^ 



0.54"' 



< 0.17(0.00) 
62+0-12 

^■"^-0.13 
14+0-34 
^•-^^-0.13 

SO+0-07 
41+0-02 

^• + 1-0.04 

44+0.07 

^•^^-0.06 

n 12+0-12 



o.83t°o;?^ 

2.74+0;i? 

3.05 
-2.62 
+ 1.48 
-0.15 
+0.24 
-0.29 
+0.10 
-0.21 
0.25 
0.11 



-2.50 
2.88 
2.18 
0.95 
0.80 



1.387 (39) 
1.150 (62) 
1.083 (50) 
1.682 (61) 
1.674 (104) 
1.528 (77) 
1.335 (64) 



Note. — Errors are with 1 a confidence. The 1 a hmits (with the best-fit values in the parentheses) are presented where the best-fit parameters are 

unconstrained. I 

h 

Regions as marked in Fig. 1. 
''Ionization timescale for the non-equilibrium ionization plasma. 
^Emission measure J neni/dV (dz,/5.2kpc)^ {cLl is the distance toward RCW89.) 

'^F^Yi ^PL ^I'c the unabsorbed energy flux in 0.5 — 5.0 keV for the thermal component and for the power-law function, respectively. They are expressed 
in a unit of ergs cm~2 g-i. 
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Table 3. Estimated plasma parameters assuming di — 5.2 kpc 



Region ^ 


7 7 b 

la X /ft 


T /" C 

V 


He 


-^thermal 


''plasma 






\cm ) 


\J ) 




ergs; 


( fl/2 \ 


1 


2.0 X 1.7 














1.5 X 1.0 


1.5 X 10^2 


510lI°o 




X 10^^ 


> 59(150) 


2 


8.0 X 2.5 














2.0 X 2.0 


1.2 X 10^^ 




9 1 +0.4 
^■-•--0.3 


X 10^6 


> 220(660) 


3 


2.2 X 2.0 














2.7 X 2.7 


5.7 X 10^2 




9 c,+1.5 


X 10^6 


10+^0 


4 


5.0 X 2.4 


5.8 X 10^2 


160^^ 


1.81^:? 


X 10^6 


22;^ 


5 


15.7 X 4.9 


7.5 X 10^3 


lOOl^ 




X 10^7 




6 


3.4 X 2.0 
11.3 X 2.7 














1.7 X 1.7 


2.0 X 10^3 


761^ 


q O+0.2 

-^■0-0.3 


X 10^^ 




7 


3.8 X 2.0 


3.0 X 10^2 


130li^ 


1 q+0.3 
^■-^-0.3 


X 10^6 


5li 



Note. — Uncertainties are all at 1 o" confidence. The 1 a limits (with the best 
values in the parentheses) are presented where the best-fit parameters are uncon- 
strained. There still remains uncertainty /, which is the filling factor or the uncer- 
tainty of the emitting volume. 

^Regions as marked in Fig. 1 

^Size of the plasma cloud. The axis length (2Za, 21^) are defined as the FWHM of 
the surface brightness. For the region 1, 2, 3, and 6, we measured the size of each 
distinctive cloud. 

^Total volume of the plasma clouds in the region. 

"^Internal thermal energy contained in the knots. 

'^Plasma age obtained from the ionization parameter and electron number density 
{net/rie). 
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RA(J2000) 



Fig. 1. — (left) An overall X-ray image of RCW89 and PSR B1509-58 in the energy range 0.4—8.0 
keV. The pulsar position is marked with a plus sign, (right) A zoom up of the interacting region 
in RCW89 The ellipses and rectangular regions were selected for succeeding spectral analysis. 
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Fig. 2. — Left: Power law fit to the jet region(labeled in Fig. 1) — Right: Model fit to the 
brightest region in RCW89 labeled "A" in Fig. 1. We used a non-equilibrium ionization model, 
"VNEI" in XSPEC, for the softer thermal component (dashed line) and a power law model for 
the hard component (doted line). 
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Fig. 3. — Spatial distribution of the temperature (Left) and the ionization time scale Uet cm~^s 
(Right). The plasma temperature increases along the "horse-shoe" in clockwise direction, while 
in contrast, Uet decreases. Uncertainties are at 1 o" confidence. 




2 3 4 5 6 7 



Fig. 4. — Left — The observed redshift for the each plasma cloud. Right — The schematic image 
of the geometry of the pulsar and RCW89. The errors are 1 a. 



